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SUMMARY 


A controller design approach for large space structures is presented, 
which consists of a primary attitude controller and a secondary or damping- 
enhancement controller. The secondary controller, which uses several 
Annular Momentum Control Devices (AMCD's), is shown to make the closed-loop 
system asymptotically stable under relatively simple conditions. The primary 
controller using torque actuators (or AMCD's) and colocated attitude and 
rate sensors is shown to be stable. It is shown that the same AMCD's can be 
used for simultaneous actuation of primary and secondary controllers. 

Numerical results are obtained for a large, thin, completely free plate model. 

INTRODUCTION 

Future utilization of space is expected to require large space structures (LSS) 
in low Earth and/or geosynchronous orbits. Examples of such future missions 
include: electronic mail system. Earth observation systems., solar power 

satellites, and space manufacturing facilities, requiring large antennas, 
antenna platforms, and solar arrays. These missions will be feasible because 
of the availability of the space Shuttle for relatively inexpensive transportation 
into low-Earth orbits. Shuttle capability can be expanded by augmenting 
a low-Earth to geosynchronous orbit transportation system. 

To establish these structures in space at minimum cost will require that 
their weight be minimized. As a result, these structures will tend to nave 
extremely low-frequency, lightly damped structural modes which are closely- 
spaced in the frequency domain. Structural parameters (i.e., frequencies. 
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damping ratios and mode shapes) are usually difficult to determine a priori. 

For these reasons, control systems design for LSS is a complex and challenging 
problem. Two types of control systems will be required for LSS: (i) large- 

angle maneuvering in order to reorient the LSS, and (ii) pointing the LSS in 
space with the required precision in attitude and shape. The objective of 
this report is to develop and investigate a controller design methodology for 
pointing control of LSS. 

The basic problems in pointing control of LSS have been well known for 
several years in the context of control of conventional spacecraft, which are 
relatively rigid, but which have sufficient flexibility to necessitate con- 
sideration in the design process. However, structural flexibility is the most 
dominant characteristic of LSS, making LSS a new class of spacecraft. Because 
of pointing requirements it is necessary to have LSS closed-loop rigid-body 
bandv/idth higher than a number of structural mode frequenices. Because of the 
high order of the LSS state vector (resulting from a large number of dominant structural 
modes), a practical controller can be designed to actively control only a 
few of the structural modes. Stability of the system is not assured with 
low-order controllers, because of control "spillover" (i.e., unwanted forcing 
of the uncontrolled or "residual" modes by the control input) and observation 
"spillover" (i.e., unwanted contribution of the residual modes to sensor outputs) 

(ref. 1). These problems were considered in references 2, 3 and several 
methods for designing reduced-order controllers based on Linear-Quadra tic- 
Gaussian (LQG) control theory were proposed and discussed in reference 3. 

Stability of the closed-loop system designed using these methods is heavily 
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dependent on the inherent damping present in the residual modes (ref. 2). 

As mentioned previously, inherent damping ratios of residual modes are 
difficult— if not impossible— to predict. Therefore, it is highly desirable 
to enhance modal damping of LSS v/ith a robust "secondary" or "damping 
enhancement" controller. The controller design methodology proposed and 
investigated in this report takes this approach. The proposed design con- 
sists of a two-level control system which includes a primary attitude controller 
and a robust secondary controller. 

Direct velocity feedback controllers (ref. 4) have been proposed in the 
literature for damping enhancement in LSS. In particular, "member damper" 
controllers (ref. 5) and low-authority structural controllers (ref. 6) with 
guaranteed Lyapunov stability, have been proposed. In this report, the use 
of several Annular Momentum Control Devices (AMCD's) for damping enhancement 
is proposed and investigated. (See ref. 7 for a description of AMCD. ) The 
secondary controller makes the closed-loop system asymptotically stable under 
certain relatively simple conditions. Primary attitude controller design is 
considered using either torque actuators or AMCD's. When the torque actuators 
are colocated with attitude and rate sensors, the closed-loop system is 
stable. It is proved that the same AMCD's can be used to accomplish 
simultaneous primary and secondary controller actuation. The overall control 
system is shown to be stable for the case when actuators are colocated with 
sensors. Numerical results are obtained for a large, thin, completely free, 
flat plate in order to demonstrate and evaluate the controller design methods. 
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MODEL INVESTIGATION 


Mathematical Models of LSS 

A large space structure is a highly oscillatory distributed parameter 
system. A class of large space structures is described by the partial dif- 
ferential equation: 



+ A. 


d 

at 


uU,t) 


+ m(4) 


^ A' lo(-6 }t) - f('4 jt) 

ar 


(1) 


where A-| is a linear operator consisting of partial derivatives of the 
deflection function u(4,t) with respect to space variables a, is 
a linear operator describing the inherent structural damping in the LSS, 
and functions m(4) and f(4,t) denote the mass distribution and the 
applied (generalized) force distribution. Assuming zero damping (A 2 = null 
operator), it is possible to obtain the following normal -coordinate representa- 
tion using the property of orthogonality of eigenfunctions of A-j, and 
using appropriate boundary conditions: 


2 m 

ij,-+ 0 .,. q, = 


(i— 1 ,2 , . . . ,°°) 


( 2 ) 


where q. is the modal amplitude of the ith mode, Uj^, k=l,2,...m represent 
the input which is assumed to consist of generalized point-forces (i.e., 
forces and moments), denotes the value of the ith "mode shape" at the 

location of kth actuator. The model consists of an infinite number of modes. 
Since real-life LSS will have some structural damping, a more realistic 
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representation can be obtained by adding the term 2p^.w^. to the left-hand 
side of equation (2), where p. and denote damping ratio and natural 
frequency of the ith mode. The "observations" or sensor outputs, which 
consist of position and angular displacements, are given by 

oo 

'^ki^i (3) 

where denotes the ith mode-shape evaluated at the location of the kth 
sensor. The modes in equation (2) also consist of rigid-body modes. It is 
customary to separate the rigid-body modes from the structural modes, and 
also to truncate the model at structural modes. For practical LSS, 
it is usually not possible to analytically obtain the model of the form of 
equations (1) and (2). The standard technique is to use the finite element 
method to obtain a model. Established computer programs such as SPAR and 
NASTRAN are available for this purpose. Finite element structural models 
also have the same form as the normal -coordinate model (i.e., eqs. (2) and (3)). 

For the purpose of this investigation, it was necessary to choose an 
appropriate model of a LSS as a medium for control law development. After 
consultation with the NASA technical monitor, a finite element model of a 
30.48 m X 30.48 m x 2.54 mm (100 ft x loo ft x o.l in.), completely free 
aluminum plate was selected for this purpose. This model was developed in 
reference 8 using the SPAR program. A finite element mesh of 24 x 24 equal 
square plate elements was used to obtain modal frequencies and mode shapes 
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(with respect to force and torque inputs)which were computed(ref. 8) for the 
first 44 modes. Table I gives the rigid-body parameters and the first 44 modal 
frequencies for this model. Values of the mode-shapes at all 625 nodes are 
given in reference 8. Rotation about only two axes (x and y axes in 
fig. 1), and translation in the z-direction were considered in the model since 
they suffice to demonstrate the .principles. Thus the LSS model is given by: 


(i) Rigid-body motion: 


Ht 

Z T. + Z R,- X f. (4) 

j=l J i=l ^ ’ 


Of 

m z = Z f . (5) 

s s ^.^1 1 

where a = (4> denotes the rigid-body attitude vector about x and y axes; 

s s s 

z denotes z-axis translation of the LSS center of mass; K and m^ denote the 

S 5 5 

two-dimensional LSS inertia matrix and LSS mass; f^. and Tj denote applied 
forces and torques (i=l ,2,. . ; J=1 ,2,. . ,ny) ; R^. denotes the coordinates 
of point of application of force f.. It should be noted that the rigid-body trans- 
lation Zg is not of interest in investigations of the LSS attitude control problem. 
Equation (5) is included here for completeness, and will be used only in equation (7). 
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(ii) Flexible motion (assuming no inherent damping) 

q + Aq = (6) 

where q is the riq-dimensional modal amplitude vector; A is a positive 
definite pq x nq matrix (usually diagonal, with entires = <l>^ and 

are n^ x Oq and nj ^ Pq mode shape matrices corresponding to force 
and torque inputs, and f and T denote vectors consisting of f^. and 
T^. , respectively. 

(iii) Sensor outputs: 

Sensed z-axis translation is given by (ignoring noise) 

^m = ^s + Rz ^ “s •*' '*’f^ (7) 

where is the coordinate (vector) of the sensor location, <|)^ is the 
Pqj^^ mode shape vector at the sensor location. Sensed LSS attitude 
(4>m>6m^ is given by 



where are the nq^i mode shape vectors at the x and y axis 

attitude sensor locations. 
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Controllability and Observability of LSS Model 
The LSS variables to be controlled are: rigid-body attitude and 

rate (a^ and a^), and structural modal amplitudes and rates (q,q). 

Control inputs are point forces or torques, and observations (sensor outputs) 
normally consist of an attitude sensor (e.g., sun sensor or star tracker), 
and an attitude rate sensor (e.g., rate gyro). Let (i=l ,2,. . . ,n^) denote 
the ith eigenvalue of A (X. = where is the frequency of the ith 
mode). The following theorem gives the necessary and sufficient conditions of 
controllability of the system defined by equations (4) and (6), with respect 
to the torque vector T (or force vector f). 

Theorem 1 .- The system given by equations (4) and (6) is controllable 
with respect to T(f) iff (if and only if): (i) the rigid-body system 

(eq. (4)) is controllable with respect to T , and (ii) for the flexible part, for 
i=l ,2,. . . ,n^, if (a) X^. is a simple eigenvalue, '^^(j>i) t 0 for some j e [1 .n^] 

(b) if X^ has multiplicity , then rank of the n^. ^ y^- block of 

(n^ ^ y.j block) corresponding to X^ is y^-. Similar results hold for force input f. 

Outline of proof .- It has been proved in reference 9 that the system 
given by equation (6) is completely controllable if conditions (iia) and 
(iib) hold. Considering the composite system given by equations (4) and (6), 
results of reference 10 can be applied in a straightforward manner to 
complete the proof. ! 

Equation (6) does not include the damping term, which, although small, 
will always be present in practical LSS. When damping is represented, 
equation (6) becomes: 
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(7) 


q + Ddi + Aq = 

where D = ^ 0 is the damping matrix. If the system is 

controllable for D = 0, it will be controllable for D > 0, for sufficiently 
small D (because of the continuity property). However, stabilizability of 
the system is more important, especially while designing LQG regulators. 

The following theorem gives sufficient condition for stabilizability. 

Theorem 2 .- The system given by equations (4) and (6) is stabilizable 
for D 0 if it is controllable for D = 0. 

Outline of proof .- By applying the results of reference 10, it can be 
proved that the system is stabilizable if condition (i) of theorem 1 is 
satisfied, and if the system of equation (6) is stabilizable. If the 
system of equation (6) is controllable, then there exists a rate feedback 
gain such that the closed-loop system (with ^ f 0) is asymptotically stable 
(ref. 6). Thus the system of equation (6) is stabilizable for D f 0. 

This completes the proof. 

Since observability is the dual concept of controllability, theorems 1 and 2 
can be used to also investigate observability. The plate model selected for the 
LSS study has a number of modes with repeated natural frequencies (multiplicity = 2). 
Using theorem 1 it is evident that at least two actuators will be required 
for controllability (necessary condition). It is straightforward to check 
controllability (and observability) with respect to given actuator (sensor) 
locations by applying theorem 1. 
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Actuator and Sensor Model s 

The controller design approach proposed in this report uses several 
Annular Momentum Control Devices (AMCD's) for damping enhancement, and torque 
actuators or AMCD's for primary attitude control. An AMCD (fig. 2) 
consists of a rotating thin rim which is suspended in three or more non- 
contacting electromagnetic actuators and driven by a noncontacting electro- 
magnetic spin motor (ref. 7). The bandwidth of AMCD sensors and actuators is 
very high— of the order of hundreds of Hertz— therefore, AMCD sensors and 
actuators can be assumed to be without phase lag. In the AMCD model, only 
X and y axis rim rotations and z-axis rim translation need be considered. 

The basic linearized AMCD equations of motion may be found in references 11 
and 12. Since AMCD's are assumed to be relatively small (of the order of 
2 m rim diameter), the rims can be assumed to be rigid. Torque actuators and 
attitude/ rate sensors (used for primary controller) are assumed to be linear and 
to have infinite-bandwidths. This assumption is justified in case of sensors 
since sensor bandwidths are expected to be several decades higher than LSS 
closed-loop rigid-body bandwidth, and modal frequencies of interest. Effects 
of finite bandwidth and nonlinearities in torque actuators are not considered 
in this report in order to be able to obtain certain fundamental results. 

A TWO-LEVEL CONTROLLER DESIGN APPROACH 
The controller design approach proposed in this report consists of a 
primary and a secondary controller. The function of the secondary controller 
is to enhance damping in LSS structural modes without attempting to control 
rigid-body modes. The secondary controller should be robust— that is, it 
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should be stable regardless of parameter inaccuracies. A secondary controller 
consisting of several AMCD's is proposed in this report. The advantages of 
using AMCD's are: (i) precise knowledge of modal frequencies and mode shapes 

is not required, (ii) the closed-loop system is Lyapunov-stable (asymptotically 
stable under certain conditions) regardless of parameter inaccuracies and 
regardless of number of modes in the model, as will be shown later, 

(iii) associated weight penalty is small, and (iv) the secondary controller 
using AMCD's imparts gyroscopic stability to LSS. The last feature may be 
useful during initial phases of deployment, assembly or initial in-orbit 
parameter estimation, before the primary attitude controller is activated or 
even designed. 

Secondary controller design using a single AMCD was investigated in 
references 13 and 14, with the latter reference containing more detailed 
stability results. The case with several AMCD's was investigated in 
reference 15. Although the sufficient conditions obtained in reference 15 for 
Lyapunov stability are easy to satisfy, those for asymptotic stability (AS) 
are difficult to satisfy in practice. Simpler sufficient conditions for AS 
are presented in this report. 

SECONDARY CONTROLLER USING AMCD'S 
Mathematical Model of LSS/AMCD's 

It is assumed that v (>_1) AMCD's are used on an LSS of mass m^, 
inertia matrix I^ (two-dimensional), and n^ structural bending modes. 

The AMCD rims are assumed to be relatively small (=2 m) in diameter; therefore, 
they can be considered to be rigid. Only x and y axis rotations and 
z-axis translations are considered, which suffices to present the principles. 
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The location of nominal rim center position of the ith AMCD in the LSS 

coordinates is (X|^, y^.)* * r^, (^3) represent the mass, rim 

radius, transverse-axis inertia matrix, and number of actuator stations for 

the ith AMCD. C-j^. represents 2 x moment-arm matrix of the AMCD 

actuator stations, and Cg^- is a lx vector consisting of all unity 

entires, e. represents the Z axis displacement of the ith AMCD rim 

center from the corresponding point on the LSS. Let a =(<(), 6 ) 

denote the LSS attitude vector about the x and y axes, and a, . = (<j) ., 0 .) 

ai '^ai* ai' 

denote the rim attitude vector for the ith rim. The actuator stations are 
assumed to produce only axial (z-axis) forces. (Only radial rim centering is 

accomplished by radial actuator forces and is of no consequence in the 
present analysis.) 

From the basic principles of dynamics the combined AMCD/LSS equations 
of motion can be written as: 

Ax + Bx + Cx = Yf (9) 


X 





a 


s» M 


"V 


nT\T 

.q ) 


where q denotes the Uq-vector of modal amplitudes of the LSS 

f = (f{, fJ,..., F^ 

F^. = f^ 2 ***‘* ^ijj,.^ i=l,2,...,v 


( 10 ) 


( 11 ) 

(12) 
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f^.|^ being the axial force at actuator station k of the ith 
AMCD. 

A is a (n.| x n^) symmetric positive definite coefficient 
matrix (where n^ = n^ + 3v + 2): 



(3v + 2) 



^nqxnq ^^i^otes the nqXn^ identity matrix 



13 


(13) 



where c = = (y,-. - x^.). = diag. (m^ ,..., ) 

and ^ ^ 


K’}i4 = ’ ^ 


B = 


V 

S W 
1=1 


1 


w. 


V 


W, W« ' * w^, 0 0 

12 VI I 


Wn 


T"""T 


w 


0 0 


0 I 0 I D 


(14) 


where D = D >0 denotes the n„ x n„ 

- q q 


LSS damping matrix. 



( 15 ) 
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being the angular momentum of the ith rim about the Z axis. 


°(3v + 2) X (3v + 2) 




( 16 ) 


Where A = > 0 is the LSS modal frequency matrix (usually diagonal, 

with squared modal frequencies as its entries). 
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(i=l,2 v) represents the A. x n mode shape matrix for actuator 

^ ^ A 

locations of the ith AMCD, I = , and rig = n^ + 3v(= in eq. (17) 

denotes equality by definition). 


Let 6^.|^ be the axial centering error at actuator station k of 
the ith AMCD. The £ x 1 centering error vector 6 is given by 


6 = ( 6 ^^,..., 


= -T^x 


(18) 


Closed-Loop Stability 
Consider a control law of the type 

f = Kp6 + K^6 (19) 

Where Kp and are real symmetric positive definite S, ^ I proportional 
and rate gain matrices. 

It was proved in reference 15 that the closed-loop system given by 

equations (9) and (19) has at least two zero eigenvalues for all Kp and 

K . The zero eigenvalues correspond to a which represents LSS rigid-body 
X* s 

attitude. Defining: 
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(20) 


P = (ctl 






X = (p^. 





( 21 ) 


Equation (9) (excluding a^) can be expressed as: 


X = A X + B f 


where 


( 22 ) 


A 


n2><n2 


n 2><2 


n2><n 



+ 



/^(3v+2)xn2\ I 

Vn^xIvrA7 I 


-L 



(23) 


17 



denotes the mxm identity matrix) 

' mxm 



( 24 ) 


It was proved in reference 15 that the closed-loop system is Lyapunov- 
stable if Kp > 0, i 0, and is asymptotically stable if: 

V 

(i) >0, > 0, (ii) Z H. 0, (iii) £ > n +3v, (iv) rank (y) = n + 3v. 

p r 1 I H H 

However, conditions (iii) and (iv) are difficult to satisfy in practice. 

In order to obtain the least restrictive (necessary and sufficient) conditions, 
it is first proved that the control law of Equation (19) is optimal with 
respect to a certain linear quadratic performance index. 

Define the matrices 

(25) 


nKpRKpP' 


Q = 


? -1 T 

pA + p yR T 
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where A - 2 diag. (0(3v+2)»^)» p > 0 (scalar), y is given by (17). 

R = pK^’\ and S = 


Theorem 3 .- The optimal control law (provided that it exists) which 
minimizes the performance function 


nKpR 


n^x£ 


J = 




Qx + 2x^ Sf + 


f^Rf) dt 


(26) 


is given by equation (19). 

Proof.- The proof can be obtained by comparing equation (11) v/ith 
the standard solution of the LQ regulator problem (eq. (26)). 

It can be shown that the Riccati matrix is given by 


nKpn"*" + 


C !) 


0 


p = p 


(27) 


0 


A 



It has been proved in reference 15 that the matrix [nKpH^ + dia^ (0,A)] > 0. 
Therefore, P > 0. 

Stability results .- Well established stability properties of the closed- 
loop optimal LQ regulator can be used for investigating the asymptotic stability 
of the system. The AMCD/LSS system consists of the two subsystems represented 
by the following set of equations: 


— A*! x^ 

(28) 

X2 = ^ 2^2 ^2^ 

(29) 


where 


/ T T T T 

3i s a^ s 1 
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(30) 


= {q^ 




(31) 



°3vx3v 1 ®3vx2 ^3vx3v 


0 

= 

+ 

°(3v+2)x3v 1 -Ai’^B, 

1 ' 1 

. = 



L 1 J 
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where and represent the top 3v x 3v submatrices of A and B, 
and L is the 3v x 5 , matrix consisting of the top 3v rows of n. 


1 

0 

q q 

. «2 = 

1 

0 

1 

1 

> 

-D 






_ 


Since ^ 3 for all i and actuators are at distinct locations, rank 
(L) = 3V. Equation (29) can be transformed into the controllable canonical 
form: 




— 

— 


— — 

# 

"1 

= 

“11 

«12 

+ 

'P 

1 

-s« 




0 

“22 _ 


0 


where the dimension of r, is n^ (< 2 n„) and that of r« is 2n - n . 

I c — q t q c 

“ir “12’ “22 ^ appropriately dimensioned matrices. 

Theorem 4 .- The optimal control law for the problem in Theorem 1 exists, 
and the system of equation (22) with the control law of equation (19) is 
asymptotically stable, iff (if and only if) all the following conditions 
are satisfied. 
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(a) the pair (Ag, Bg) is stabilizable. 


(b) Z W. 0 
i=l ^ 


(c) 


B-|±2jn^.Ai 


rank 


an±2jfliln 


c c 


1= + 3v + 2 




for all i (i-1,2,..., v) for which / 0, where = spin velocity of 
the ith AMCD rim. 

Proof.- From linear-quadratic optimal regulator theory, it is well known 
that the control law of equation (19) makes the closed-loop system asymptotically 
stable if and only if (1) (A, B) is stabilizable, where A = A - B R”^ S^, 

and (2) (A, C) is detectable, where Q = C c"*" (ref. 16). Also, P is 

positive definite if and only if (A, C) is observable; thus statement (2) 
above is true since P > 0. Stabilizability of (A, B) is equivalent to 
that of (A, B) (ref. 17), which is equivalent to (1) stability of 022 
and (2) controllability of the composite system consisting of equation (28) 
and the equation: 
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r*! = ctn + ipf 


This composite system is controllable if and only if (ctii> is 
controllable and rank (c) = n^, + 6v + 2 where 



and where is the ith eigenvalue of (ref. 10). Assuming the AMCD 
rims to be thin it can be verified that the eigenvalues of A-j are: 

0, ±2jn. , where is the spin velocity of the ith rim. The proof can be 
completed by using elementary matrix operations. 

These necessary and sufficient conditions are not straightforward to 
apply. The following corollary gives sufficient conditions for asymptotic 
stability which are more easily applicable. 

Corollary 1 .- The system of equation (22) with the control law of 
equation (19) is asymptotically stable if all the three conditions given 
below are satisfied: (a) the LSS structural model, i.e., the pair (Ag, § 2 ) 

V 

is stabilizable, (b) Z Uj 0, and (c) the LSS does not have an undamped 

i=l 

structural mode with frequency 2J2. (i=l,2,..., v). 
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Proof.- SincG is not an GigenvaluG of <^ii» It romains to bG 

provGd only that (ThGorGin 2) 


rank ± j2fl^.A^ 


= 3v + 2 

(33) 

for i 0. 




Or equivalently that 




rank ± j2n^. | 

1 

(S 

= 3v + 2 

(34) 

By applying elementary row operations and using the fact that rank 

(Cii) = 2, (i=l,2 v) , the rank of the matrix in equation (34) can be shown to 

be 3v + 2 for f 0. Thus from Theorem 4, the system is asymptotically 

stable. 


ThG following corollary gonoralizos thG sufficiont conditions which wore 
dorivGd in roforoncG 15. It should bG notGd that thay ara too rastrictivG 

sincG thay raquira a larga numbar of actuators, and ara givan hara only for 
complatanass. 

Corollary 2.- Tha systam in aquation (22) with tha control law of 

aquation (19) is asymptotically stable if all four conditions given below are 
satisfied; 
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(a) The LSS structural model, i.e., the pair (Ag, § 2 ) 

V 

(b) Z W^. f 0; (c) 5,^ n + 3v, and (d) rank (y) = n + 3v. 

1 ^ 

Pt^oof -- Since rank (y) = n. + 3v, rank 
of Theorem 4 is satisfied if for X = ±j2Q^., 


L 


= n^ + 3v. 


is stabilizable; 


Condition (c) 


rank 


± XA^ 



= n^ + 3v- + 2 



= 3v, condition (c) of Theorem 4 is satisfied 
B^-XA^ are independent for X = + 2jJ2^.. 


Since rank 

if the two top rows o1 
It can be verified that this is indeed the case. 

The sufficient conditions in corollary 1 are less restrictive than those 
in corollary 2 because the former do not require a large number of actuators. 

It should be noted that, in addition to strutural damping enhancement, 
the secondary controller also imparts gyroscopic stability to the LSS 
(rigid-body attitude). This feature should be useful during assembly and 
during initial operations (such as parameter identification) in a newly 


deployed or assembled LSS, before the primary controller is activated, or 
even designed. 
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PRIMARY AHITUDE CONTROL SYSTEM 

The secondary control system increases modal damping and thus controls 
the shape of LSS. It also aids in primary controller design by reducing the 
effect of "spillovers." In this report, primary attitude controller design 
is considered using torque actuators and AMCD's. A primary attitude control 
system using torque actuators is considered first. 

Primary Controller Using Torque Actuators 
Assuming that the primary attitude control is accomplished using 
a (^1) two-axis torque actuators distributed on the LSS, the LSS equations of 
motion (without secondary controller) are given by: 


KK 

s s 


B3X5 


+ C X 
s s 


a 

= E 
i=l 


Yj. "T. 
'ti 1 


(35) 


where x^ = (a^^, q^)^, A^ 
Cg = diag. (0, A). 


diag. (Ig, I^ ), = diag. (0, D), 

q q 



CVJ 
X 
OJ 
1— 1 

1 


1 

•r— 

X 

1- 

1 

Yti = 


• T,- 





T . 


rrn 


yi 


(36) 
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being the 2 x "torque. mode-shape matrix" corresponding to 
location of actuator i. and are the x and y axis torques. 

Colocated actuators and sensors .- Assuming that a, two-axis attitude 
and rate sensors are also located at the same points as the torque actuators, 
the measured attitude vector and the measured attitude rate vector 

“mri location i (ignoring noise) are given by: 


a . 


mi 



(37) 


a 


mri 



(38) 


Denoting r = ,. . . , 
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“mr " l^°‘mrl 
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equations (35), (37), and (38) can be written as; 
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Choosing the control law 


7 = _(G a + G a ) (421 

' p m r mr^ ' 

where 6p and are 2a x 2a positive definite symmetric matrices. 

For example, a simple control law for controlling the rigid-body attitude 
would be: 


)<4i 

(43) 

for i=l,2,..., a, where and (k = x,y) denote the desired rigid- 

body damping ratio and natural frequency. In this case Gp and G^ will 
be block-diagonal matrices. 

The closed-loop equations become: 





( 44 ) 
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Finally, the AMCD loops may be closed using equations (18) and (19). 

Theorem 5 .- The system of equation (44) is stable in the sense of 
Lyapunov if Gp > 0, ^ 0. 

Proof .- Considering a Lyapunov function: 



(since B = B + TG > 0). 

o o I 

Thus, the primary control law using colocated torque actuators and 
attitude and rate sensors gives a stable closed-loop system. The controller 
is robust because it is stable regardless of parameter inaccuracies. In 
addition to controlling rigid-body modes, the "colocated" primary controller 
also increases damping in some of the structural modes, depending on the 
locations of torque actuators. The secondary controller consisting of 
AMCD's can be added in order to enhance modal damping and improve the overall 
performance. When the secondary controller is included, it can be shown that 
the overall closed-loop system is Lyapunov stable if Gp >0, G^ ^ 0, 

K > 0 and K > 0. 
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Noncolocated actuators and sensors.- Stability of the closed-loop 


system is no longer guaranteed when actuators and sensors of the primary 
controller are not colocated. This is because of control and observation 
"spillovers." In reference 3, the following controller design approaches 
based on LQG control theory were investigated: 

1. Truncation: In this method, the residual (uncontrolled) structural 

modes are merely ignored in the design process. 

2. Modified Truncation (or Model Error Sensitivity Suppression, 

Reference 18): The effect of control input on selected residual modes is 

« 

included in the performance function in a quasi-static sense. 

3. Use of Higher-Order Estimator: In this method, the state estimator 

estimates more modal amplitudes and rates than are fed back. 

4. Selective Modal Suppression: Observation spillover for selected 

residual modes is reduced or eliminated by using such devices as phase- 
lock-loops. 

5. Polynomial Estimators: Observation spillover is explicitly 

estimated by representing it as a polynomial in time. 

These methods typically result in a controller which consists of a linear 
regulator and a Kalman-Bucy filter for state estimation. It was reported 
in reference 3 that, when a sufficient number of sensor outputs are available, 
better results are obtained using direct sensor feedback (DSF) than using 
state estimators. Of the methods considered in reference 3, method 5 was 
reported to be unsatisfactory. Preliminary numerical results obtained by 
applying methods 1 to 4 to the plate mode (discussed previously) indicated 
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that the modified truncation { or Model Error Sensitivity Suppression (MESS)} 
method was the most promising method. Therefore, only the truncation method 
(being the most straightforward) and MESS method (being the most promising) 
are discussed in this report. As stated earlier, DSF was reported to result 
in better performance; therefore DSF is used in this report instead of 
state estimators. As a simple but important special case, it is assumed in 
the following analysis that the primary controller controls only rigid-body 
attitude, without attempting to actively control any structural modes. 

This assumption is justified to some extent because structural modes will be 
controlled by the secondary controller. In order to simplify implementation, 
the primary control law is constrained to require feedback of only measured 
attitude and rate (measured attitude and rate includes contributions of 
rigid-body and structural modes). It should be noted that this special case 
is included only for the purpose of demonstration and that primary controller 
design using truncation or MESS methods is not restricted to active control 
of rigid-body modes only. 

MESS method - a special case: Assuming that the primary controller 
is to be designed to control only rigid-body modes [=(4>5»05)]» and 
for the case with a single 2 -axis torque actuator, consider the | 

control law given by 
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where (i = x,y) denote the desired rigid-body bandwidth and 

damping ratio, and denote the measured attitude and rate vectors; 

= l/2u)^.. If control spillover were absent (i.e., if the actuators could 
force only rigid-body modes), this control law would minimize the performance 
function 
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(48) 


This can be verified along the lines of reference 6. However, because of 

flexibility, the control input also forces structural modes. The model error 

sensitivity suppression method requires augmentation of a quadratic function 

of the forcing terms corresponding to a few selected structural modes, 

to the performance function in equation (48). If modes k.j,...kp are 

selected for augmentation, denoting the mode shapes (for the kth mode)! at the actuator 

location by and (j>y|^, the modified performance function becomes 
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(49) 


J = 



'^xk'^’yk 

^yk2 


T 


where Qj^ is a positive weighting coefficient. The effect of addition of 
this term to the performance function is to modify the control weighting 
matrix. The optimal control law for this type of performance index requires 
feedback only of measured attitude and rate. A modified state estimator can 
be similarly designed for generating estimates of rigid-body attitude and 
rate. 


Primary Attitude Control Using AMCD's 
Instead of torque actuators, AMCD's can be used for primary 

attitude control. In particular, the AMCD's used for implementing a secondary 

controller can be simultaneously used for the actuation of the primary attitude 

controller. Primary attitude control is accomplished by torquing against 

AMCD momenta. In this dual control mode, however, the relative rotation 

angles between the LSS and the AMCD's (i.e., - a.s> i=l»2,.., v) |cannot be 

controlled simultaneously with LSS rigid-body attitude (controllability 

of AMCD/rigid-body modes is discussed in reference 12). This entails that 

the control law of equation (19) for the secondary controller should be 

redesigned to exclude feedback of (a - aj. However, the rates (a - a ) 

s s' 

must be zero in steady-state and must be fed back. That is, the position 
feedback gain must be redesigned to control only , the motion of the AMCD 
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rim centers relative to the LSS, and AMCD rim transverse rotation angles are 

allowed to be nonzero in steady-state. Provided that relative angles, 

«=> - are sufficiently small (so that the actuator gap limits are not 
®i ^ 

exceeded), the AMCD's can be used in this dual control mode. 

Columns of the 2 x matrix C-j^. are given by (y^*j» - 

j=l,2 Since actuators of each AMCD are located along a circle, 

columns 1 and 2 are linearly independent, and columns 3 through can 
be expressed as linear combinations of the first two columns. That is, 

C-j^. can be expressed as: 


I ^i^ 

where is the 2x2 matrix formed by the first two columns of 
C-|^, I 2 is the 2x2 identity matrix, and r^. is a 2 X(£^. - 2) 
matrix. If Kp is designed as follows: 



where Kp is a (5--2v)x(Jl-2v)j positive definite matrix, and diag. ( ) denotes 

a block-diagonal matrix, it can be verified that the resulting coefficient matrix 

multiplying (a^ - aL..., ~ al)^ in equation 19 (after substituting for <5 

a-| s % ^ 

from equation 18) is zero. Defining 
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the resulting closed-loop system is the form: 
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where a-j , 0-2 are appropriately dimensioned matrices. 

Theorem 6 .- The system defined by equation (53) is stable in the sense 
of Lyapunov if Kp > 0 and ^ 0. 

Outline of proof .- The proof is similar to that of Theorem 2 (part a) 
in reference 15, except that it is additionally necessary to prove that 
^2*^p^2^ > 0 (positive definite) for the specially structured Kp of 
equation (51). This can be proved by establishing that z^C 2 KpC 2 ^z can 
be zero for some z ^ 0 if and only if 


[Fi''', -I] cJi = 0 (i=l,2,..., v) 


(54) 


Using the fact that actuators for each AMCD are located along a circle, it 
can be proved that equation (54) cannot hold. Therefore, 



This type of secondary controller essentially provides rate feedback 
for modal damping enhancement. The only function of position gain is to 
keep the AMCD rim centers at their nominal positions. Additional force 
commands can now be superimposed on the electromagnetic actuators in order to 
produce the desired primary control torque for controlling a^. Since AMCD's 
are small compared to the LSS, the effect of control moments generated in this 
manner would approximate point-torque actuators. If an attitude sensor and 
a rate sensor are located on the LSS at the nominal position of the center of 
each AMCD, this configuration would approximate colocated torque actuators 
and attitude/rate sensors, and should therefore have the associated Lyapunov- 
stability property. In this configuration, the AMCD's must have sufficiently 
large momenta in order to exert the magnitude of torque required to achieve 
the desired rigid-body bandwidth (without exceeding the electromagnetic 
actuator gap limits). Separate AMCD's may also be used for primary control 
actuation. The position gains for the rim suspension control system should be 
structured as discussed above to retain the closed-loop stability properties 
of the structural modes. For orbital applications it will be necessary to 
gimbal the AMCD's for primary controller actuation. 

NUMERICAL RESULTS 

For the purpose of demonstration of the primary and secondary controller 
design methods, the 44-mode finite element model of a 30.48 m x 30.48 m x 2.54 mm 
(100 ft X 100 ft X 0.1 in.), completely free, aluminum plate (discussed earlier) 
was used. The inherent damping ratios of all the structural modes were 
assumed to be zero. The AMCD's were chosen to have rims having 1.79 m 
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diameter and 34 kg ma^S* suspended in four equally spaced electromagnetic 

actuators and spinning at 5QQQ RPfl. Secondary controller design was considered 
first. 

Three AMCD's, centered at coordinates (4.44 ra, -8.26 m) , (-12.06 m, 

-0.635 m). and (14.6 m, -14.6 ra) were used for the secondary controller (in 
the coordinate system with axes parallel to the plate edges and origin at the 
plate center). Gain matrices Kp and were assumed to be diagonal, with 
entries kp , corresponding to the ith AMCD. Keeping the position gains 

constant at kp^ = 146 N/m, kp 2 = 14.6 N/m, kp 3 = 14.6 N/m, rate gains (k^.) 
were increased for the three AMCD's, starting with zero rate gains. With 

*^r2 ~ *^r3 ~ *^rl fi'f'st increased starting from zero. Keeping 

k 1 constant at its nominal value of 5636 N-sec/m, and with k =0 k 

r3 ’ r2 

was next increased from zero. Finally, with k^^ and k ^2 ^i><ed at their 
nominal values (k^^ as above and k ^2 = 2050 N-sec/ra) , k ^3 was increased from 

0. As shown in the root loci of figure 3, damping ratios of the structural 
modes increase, the lowest damping ratio being 0.07 for nominal gains (i.e., 
position gains, k^-j , k ^2 above, and k ^3 = 7174 N-sec/m). The root loci 
turn back towards the imaginary axis for higher rate gains. Although only the 
first seven modes are shown in figure 3, all modes exhibit similar behavior. 
Addition of each AMCD generally improves the closed loop damping ratios. Damp- 
ing enhancement in different modes depends on the values of the mode shapes at 
the AMCD actuator locations. 

Primary controller design using two-axis torque actuators and attitude/rate 
sensors was next considered. As discussed in the previous section, the primary 
controller was to be designed to control only rigid-body attitude, without 
attempting to actively control any structural modes. 
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Secondary controller gains (when used) were set at their nominal values 
given above. Numerical results were obtained using either a single torque 
actuator (2-axis) located at the center of the plate, or three torque actua- 
tors placed at the locations of AMCD centers given previously. Both colocated 
and noncolocated actuators/sensors cases were considered. For the noncolocated 
case, a single torque actuator located at the plate center was used, and the 
attitude and rate sensors were located at (15.24 m, 0 m). The "evaluation" or 
"truth model" for the noncolocated case was assumed to consist of rigid-body 
modes and the first seven structural modes. Numerical results were obtained for 
the following cases. 

1. single primary actuator with noncolocated sensors - truncation method; 

2. single primary actuator with noncolocated sensors - MESS method; 

3. single primary actuator with colocated sensors; and 

4. three distributed primary actuators with colocated sensors. 

Numerical results for cases 1 through 4 were first obtained without the 

secondary controller, and then with the secondary controller (Kp and being 
at their nominal values). The objective was to vary the primary controller 
position and rate feedback gains (or weighting coefficients in the case of 

MESS method) in order to get the highest rigid-body closed-loop bandwidth 
with the restriction that the rigid-body damping ratio (p^jj) does not fall below 
0.5 . Closed-loop damping ratios of the structural modes must also be reasonably 
high in order to obtain satisfactory shape/vibration control. Figure 4 shows a 

bar graph of the maximum rigid-body bandwidth achieved such that 

P ^5 >0.5, and of which represents the lowest closed-loop damping ratio 

for structural modes. Since the purpose of these computations was to gain some 
insight into performance of the methods discussed, formal numerical optimization 
routines were not used for obtaining rather, it was accomplished by 

varying the parameters mentioned above and observing the trends. When the 
secondary controller was not used, computations revealed that it was 
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not possible to obtain a stable design for the noncolocated case with either 
truncation or MESS methods, as indicated by zero value of w in figure 4. 
With the secondary controller included, however, stable designs were obtained 
for the noncolocated case, vn'th the MESS method causing slight improvement 
over the truncation method. For the colocated cases, although closed-loop 
stability is guaranteed, it is possible to have zero closed-loop damping for 
some structural modes (with no secondary controller), as indicated by the 
results for case 3. In addition, because of interaction of structural modes, 
the maximum achievable rigid-body closed-loop bandwidth is also limited with 
this type of control law, although all eigenvalues are guaranteed to be in the 
closed left-half of the complex plane. The highest w (about 0.05 rad/sec) 
was obtained using three distributed torque actuators with colocated sensors, 
when the secondary controller was included. It should be possible to 
increase it further by using additional AMCD's in the secondary controller. 

In all the cases considered, addition of secondary controller caused 
significant improvement. For investigating the use of the same AMCD's in 
primary and secondary controllers, preliminary numerical results were obtained 
for the case where attitude and rate sensors were located on the LSS at the AMCD 
center nominal locations. There was very little difference in the damping- 
enhancement characteristics in spite of using the specially structured K 

P 

matrix. The overall closed-loop system was asymptotically stable even though the 
actuators/sensors in this case are only approximately colocated. However, it 
will be necessary to have larger total angular momentum in order to get a 
rigid-body closed-loop bandwidth of the order of 0.05 rad/sec. 
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PERFORMANCE EVALUATION 


After primary and secondary controllers are designed to obtain satisfactory 
closed-loop dynamics (based on the known model parameters), the next step is 
evaluate the closed-loop performance in the presence of disturbances (such 
as gravity gradient, geomagnetic torques, solar pressure, etc.) and sensor/ 
actuator noise. Gravity gradient and geomagnetic torque are slowly varying 
disturbances which need not be considered for investigating dynamic performance. 
(Depending on the LSS orbital configuration, they must be compensated for 
by using slowly varying input bias torques; however, they were not considered 
in this report since the scope of this investigation is limited to fine- 
pointing control over relatively short segments of time.) In the presence of 
disturbances and sensor/actuator noise, which can be described by zero-mean 
white noise, the overall closed-loop equations can be expressed as: 

* V (55) 

where x^ is the n-dimensional state vector and v is a s-dimensional zero- 
mean white-noise with covariance intensity matrix V, which represents 
disturbances as well as sensor/actuator noise. A^ is the closed-loop 
system matrix, and is the effective noise input matrix. Bandlimited 
white noise can also be represented by this formulation by incorporating the 
associated filter dynamics in A and B . The covariance of x. evolves 
according to the equation: 
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(56) 


i = A^Z + E aJ + V bJ 

Where E =£[x^xj]. 

In stea(dy-state, Z = 0, and the resulting Lyapunov matrix equation can 
be solved for E (steady-state value of E) using one of the many available 
numerical methods. The method used in this report is that given in reference 19 
since it was found to have good convergence properties. Closed-loop performance 
can be evaluated by examining elements of E. The x and y axis RMS 
pointing error at a particular point the LSS surface is given by: 

“x = (57) 

where the total attitude angles (<J> about x-axis and 0 about y-axis) 
are given by cf) = and 0 = The RMS (1 a) and 3a-errors at 

various points of interest on the LSS can be computed in this manner. It should 
be noted that it is necessary to have knowledge of the LSS parameters s 
disturbances and sensor/actuator noise in order to obtain reliable error 
estimates in this manner. 


CONCLUDING REMARKS 

A controller design approach for large space structures was presented, 
which consists of a primary attitude controller and a secondary or damping- 
enhancement controller. The primary controller uses either torque actuators 
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or Annular Momentum Control Devices (AMCD's) to control rigid body modes 
(and possibly some structural modes). The secondary controller uses several 
AMCD's and is shown to make the closed-loop system asymptotically stable 
under relatively simple conditions, regardless of parameter inaccuracies 
and number of structural modes in the model. The primary controller using 
torque actuators and colocated attitude and rate sensors is 
stable in the sense of Lyapunov (using positive definite feedback of measured 
attitude and nonnegative definite feedback of measured rate). A method for 
structuring the position feedback matrix was given, which permits the use of 
the same AMCD's for the actuation of primary and secondary controllers. 
Generic stability results, as well as numberical results obtained for a 
large, thin, completely free plate indicate that a control system consisting 
of a primary controller using several colocated actuators and sensors 
distributed on the LSS, and a secondary controller using several AMCD's, 
holds significant promise. 
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TABLE I.- COMPUTED NATURAL FREQUENCIES OF 
PLATE STRUCTURAL MODEL 


MODE 

' FREQ(RAD/SEC) 

FREQ(HZ) 

1 

•54999E-01 

.875348-02 

2 

.80024E-01 

.127368-01 

3 

.99111E-01 

♦157748-01 

A 

.•14211E+00 

.226188-01 

5 

•14211E+00 

.226188-01 

6 

•24948E+00 

.397078-01 

7 

•24948E+00 

.397078-01 

6 

,260C8E+00 

.413928-01 

9 

•28286E+00 

.450188-01 

10 

.31515E+00 

.501578-01 

11 

.43068^+00 

.685458-01 

12 

.430688^-00 

.685458-01 

13 

•47824E+00 

.761148-01 

14 

.50003E+00 

.795838-01 

15 

.53689E+00 

.854498-01 

16 

•536B9E+00 

.854498-01 

17 

.62422E+00 

.993478-01 

16 

.65958E+00 

.104978+00 

19 

.688088^00 

.109518+00 

20 

•80973E+00 

.128878+00 

21 

.80973E+00 

.128878+00 

22 

•83371E+00 

.132698+00 

23 

. .87377E+00 

. 139068+00 

24 

•87982E+00 

.140038+00 

25 

.879B2EfOO 

.140038+00 

26 

.99216E+00 

.157918+00 

27 

.992168+00 

.157918+00 

28 

.114838+01 

.182758+00 

29 

•119228+01 

.189748+00 

30 

.119968+01 

.190938+00 

31 

.121948+01 

.194078+00 

32 

.122508+01 

.194978+00 

33 

.125328+01 

.199468+00 

34 

.125328+01 

.199468+00 

35 

.137428+01 

.218728+00 

36 

.140828+01 

.224128+00 

37 

.14871E+01 

. 236688+00 

38 

.148718+01 

.236688+00 

39 

.160598+01 

.255588+00 

40 

.160598+01 

.255588+00 

41 

.169728+01 

.270128+00 

42 

.169728+01 

. 270128+00 

43 

.171118+01 

.272338+00 

<»4 

.175238+01 

.278895+00 
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Figure 1.- Location of joints for plate structural model 
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IMAGINARY 



Without Secondary Controller 



CASE NUMBER 

Highest Closed-Loop Rigid-Body Frequency 



CASE NUMBER 

Lowest Closed-Loop Structural Damping Ratio 


Case 1. Non-colocated primary actuator (single)/sensor - 
truncation method 

Case 2. As above, but modified truncation method 

Case 3. Colocated primary actuator (single)/sensor 

Case 4. Colocated primary actuators (threel/sensors 

Figure 4.- Numerical results for primary attitude controller 
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